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1 Introduction 



The recently reported measurements of time dependent CP asymmetries in 
B ^Ks decays * by BaBar [3] 

Sm{2(3{^Ks))BaBar = -0.19l°j^ ^ q qq 

and Belle [4] 

sm{2P{^Ks)) Belle = -0.73 ± 0.64 ± 0.18 (2) 

result in the error weighted average 

sin(2/3($i^5))a.e = -0.39 ± 0.41 (3) 

with errors added in quadrature. The value in (2) corresponds to the coeffi- 
cient of the sine term in the time dependent CP asymmetry [6], see Section 
IV. Belle also quotes a value for the direct CP asymmetry Acp = —C^k, i.e., 
the cosine term, C^k — —0.19 ± 0.30 [4, 5]. Although there are at present 
large theoretical uncertainties in calculating strong phases, we still examine 
direct CP asymmetry in the paper in order to obtain qualitatively feeling for 
effects of new physics on CP violation. 

In the SM the above asymmetry is related to that in B J/'^Ks [7]- [10] 

by 

s\n{2(3{^K)) = sin(2/3(J/*X))+0(A^ (4) 

where A ~ 0.2 appears in Wolfenstein's parameterization of the CKM matrix 
and s\n{2l3{J/-^Ks,L))worid-ave = +0.734 ±0.054. Therefore, (3) violates the 

tThe 2003 new results are: % = -0.96±0.50tHi by Belle[l] and +0.45 ±0.43 ±0.07 
by BaBar [2]. 
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SM at the 2.7 a deviation. The impact of these experimental results on the 
validity of CKM and SM is currently statistics hmited. Future prospects at 
the S-factories are that the statistical error a^Ks{st<^t) can be expected to 
improve roughly by a factor of three with an increase of integrated luminosity 
from O.lab^^ to lab^^ [11] and it will take some time before we know the 
deviation with sufficient precision to draw final conclusions. 

However, the possibility of a would-be measurement of sm{2P{^Ks)) — 
—0.39 or a similar value which departs drastically from the SM expectation 
of (4) has attracted much interest to search for new physics, in particular, 
supersymmetry, two Higgs doublet model (2HDM), and model-independent 
way [12, 13]. In the paper we consider the decay B ^Ks in a model III 
2HDM. It is well-known that in the model III 2HDM the couplings involv- 
ing Higgs bosons and fermions can have complex phases, which can induce 
CP violation effects, even in the simplest case in which all tree-level FCNC 
couplings are negligible. The effect of the color dipole operator on the phase 
from the decay amplitudes, A$ = argf(^), in b ^ sss in the model III 
2HDM has been studied in the second paper of ref. [12] by Hiller and the 
result is A$ < 0.2 which is far from explaining the deviation. We would 
like to see if it is possible to explain the deviation in the model HI 2HDM 
under all known experimental constraints by extending to include the neutral 
Higgs boson (NHB) contributions and calculate hadronic matrix elements to 
the as order. Some relevant Wilson coefficients at the leading order (LO) 
in the model HI 2HDM have been given [14]. Because the hadronic matrix 
elements of relevant operators have been calculated to the q;^ order [15], we 
can obtain the amplitude of the process to the as order if we know the rele- 
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vant Wilson coefficients at the next to leading order (NLO). In the paper we 
calculate them at NLO in the model III 2HDM. Furthermore, as pointed in 
ref. [13] the NHB penguin induced operators contribute sizably to both the 
branching ratio (Br) and time dependent CP asymmetry S^k in supersym- 
metrical models. In the paper we calculate the Wilson coefficients of NHB 
penguin induced operators in the model III 2IIDM. Our results show that 
in the model III 2HDM, the CP asymmetry S^k can agree with the present 
data, S^k — —0.39 ± 0.41, within the la error. Even if the S,f,k is measured 
to a level of —0.4 ±0.1 in the future, the model III can still agree with the 
data at the 2a level. And the direct CP asymmetry can reach about —20%. 

The paper is organized as follows. In section II we describe the model III 
2HDM briefly. In section III we give the effective Hamiltonian responsible for 
B — > (j)Ks in the model. In particular, we give the Wilson coefficients at NLO 
for the operators which exist in SM and at LO for the new operators which 
are induced by NHB penguins respectively. We present the decay amplitude 
and the CP asymmetry S(f,K in i? — (pKg in Section IV. The Section V is 
devoted to numerical results. In Section VI we draw our conclusions and 
present some discussions. 



2 Model III two-Higgs-doublet model (2HDM) 



In model III 2HDM, both the doublets can couple to the up-type and down- 
type quarks, the details of the model can be found in Ref. [16]. The Yukawa 
Lagrangian relevant to our discussion in this paper is 
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sin a + cos a 
V2 



+ 



Pv^kmKI - 7^) - V^c^kmI^KI + 7^)1 U , 



-H-D 

where U represents the mass eigenstates of u, c, t quarks and D represents the 
mass eigenstates of d, s, b quarks, Vckm is the Cabibbo-Kobayashi-Maskawa 
matrix and the FCNC couphngs are contained in the matrices The 
Cheng-Sher ansatz for ^'^'■^ is [16] 



Si" = A: 



(6) 



by which the quark-mass hierarchy ensures that the FCNC within the first 
two generations are naturally suppressed by the small quark masses, while a 
larger freedom is allowed for the FCNC involving the third generations^. In 
the ansatz the residual degree of arbitrariness of the FC couplings is expressed 
through the \ij parameters which are of order one and need to be constrained 
by the available experiments. In the paper we choose ^^'^ to be diagonal 
and set the u and d quark masses to be zero for the sake of simplicity so that 
besides Higgs boson masses only Ajj, i=s, c, b, t, are the new parameters and 

will enter into the Wilson coefficients relevant to the process. 

^ Model III 2HDM has a remarkably stable FCNC suppression when one evolves the 
FCNC Yukawa coupling parameters by the RGB's to higher energies[17]. 



(5) 
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3 Effective Hamiltonian 

The effective Hamiltonian for cliarmless B decays witli AB = 1 is given by 
[13, 20] 



VZ ^ \ i=3,...,16 



Qr, + Csg Qsg + Q'r-, + C^, Q^, j + h.c. (7) 
Here Qi are quark and gluon operators and are given by 

Ql = iSaPp)v-AiPpba)v-A, Ql = iSaPa)v-AiP/3bp)v-A, 

^3(5) = {Saba)v-A^{Ql3(ll3)v-i+)A, ^4(6) = (^a^/s) V-A E(^/3^")^-(+)^' 

g g 
3 3 

^7(9) = ■7:{Saba)v-A'^eq{qi3qi3)v+(-)A, QsiW) = ■;j(Sa^/3)y-ylE^'/(^/3^a)v'+(-M' 

Qii(i3) = {sb)s+p}_^ {qq)s-i+)P: 

q 

^12(14) = \SiOj)s+p 2^ {qjqi)s-{+)p ■, 

q n^b 

gi5 = sa^^il + 75) E ^^^..(1 + 75) ? , 

Q16 = a^^il + 75) E -^Qj + 75) ?^ 

q 1Tlb 

g 

Qsg = ^mbS^a^-^G^^^il + ^,)bp, (8) 

where {V ± A){V ± A) = 7/^(1 ± 75)7^(1 ± 75), {qiq2)s±p = ± 75)^2, 
q = u,d, s, c, b, Cq is the electric charge number of q quark, Aa is the color 
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SU(3) Gell-Mann matrix, a and are color indices, and Fp,y are the 

photon [gluon] field strengths. The operators Q^s are obtained from the 
unprimed operators QjS by exchanging L R. In eq.(7) operators Qi, 
ill,.. .,16, are induced by neutral Higgs boson mediations [13]. 

The Wilson Coefficients Ci, i=l,...,10, have been calculated at LO [20, 14]. 
We calculate them at NLO in the NDR scheme and results are as follows. 



Ci(Mw 

C2(Mw 

C3(Mw 

C4(Mw 

C5(Mw 

C6(Mw 
C7(Mw 

C8(Mw 

C9(Mw 

Cio(Mw 



11 Qi,(Mw) 

2 



= 1 



11 a,(Mw) 
6 



An 



35 q; 
18 4^' 



247r 



Stt 



a 
Gtt 

0, 
a 
Gtt 
0, 



4Co(xt) + i^o(2:t) 



4Co(xt) +L>o(a;t) + 



1 



sin 9w 



[lQBo{xt) - AC^{xt)) 



A{x^) A{y) 



2 6 
D{xt) D{y) 



(9) 
(10) 

[2Bo{xt) + Co{xm 

(12) 

(13) 

(14) 
(15) 
(16) 
(17) 
(18) 
(19) 

(20) 
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where Xt — rn^/M^, and y — m^/M'^±. Here the Wilson coefficients Cy-y and 
Csg at LO which are given in ref. [14] have also been written. The Wilson 

coefficients Cj^ and Cgg at NLO in SM have been given but they at NLO 



in model III 2HDM have not been calculated yet. Because we calculate the 
decay amplitude only to the as order it is enough to know them at LO. Here 

^x2 + .5,r-7 {2,:r- -2.v)\nx- 



A{x) 
B{y) 

D{x) 

m 

Coixt) 
Do{xt) 



X 



y 



X 



12{x-lf 2(x-l)4 
5y-3 {3y-2)\ny 



L12(|/-1)2 
— 5x — 2 



+ 



6(1/ -1)3 
3a; In a; 



A{x-iy 2{x-l) 



I/-3 



+ 



1 

4 
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A{y - ly 2iy - If 
Xt , Xt In Xt 



l-Xt {Xt-l)\ 

Xt — Q 2)Xt + 2 



Inxt 
+ 



Xt - 1 (xt- 1)2 

IQxf + 25x1 x?(5a;| - 2xt - 6) 



18{xt - ly 



(21) 
(22) 
(23) 
(24) 
(25) 
(26) 

Inxt, (27) 
(28) 



and 



Eo{xt) 
Eo{xt) 



2 Xt{lS - llxt - xf) xj{15-16xt + 4:X 

■3^""^*+ 12(1 -x,)^ 



6(1 - xtY 



Eoixt) - - 



^ +— rrlny^. 



*Mnx^?9) 
(30) 
(31) 



36(1 -y)3 ' 6(1 -y)4 

The Wilson coefficients Cj, i=ll,...,16, at the leading order can be ob- 
tained from Cqi and Cq2 in Ref. [19]. The non- vanishing coefficients at niw 
are 



Cn{Mw) 



a nib 



{Cqi - Cq2) 
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(32) 



47r itItXj 

We shall omitted the contributions of the primed operators in numerical 
calculations for they are suppressed by ^ in model III 2HDM. 

For the process we are interested in this paper, the Wilson coefficients 
should run to the scale of 0(m6). Ci — Cio are expanded to 0{as) and 
NLO RGEs should be used. However for the Csg and Cj^, LO results should 
be sufficient. The details of the running of these Wilson coefficients can be 
found in Ref . [20] . The one loop anomalous dimension matrices of the NHB 
induced operators can be divided into two distangled groups [23] 



and 
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Oi5 
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-16 





1/3 


-1 


Ol4 


-6 
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-1/2 


-7/6 


Ol5 


16 


-48 


16/3 





Ol6 


-24 


-56 
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-38/3 



(33) 



(34) 



For Q'^ operators we have 

Because at present no NLO Wilson coefficients i=ll,...,16, are available 
we use the LO running of them in the paper. 
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4 The decay amplitude and CP asymmetry 
in Bl ^ 4,Ks 

We use the BBNS approach [18] to calculate the hadronic matrix elements 
of operators. In the approach the hadronic matrix element of a operator in 
the heavy quark limit can be written as 

< (t>K\Q\B >=< 4>K\Q\B >f [l + 5^r„<], (36) 

where < (f)K\Q\B >f indicates the naive factorization result. The second 
term in the square bracket indicates higher order as corrections to the matrix 
elements [18]. We calculate hadronic matrix elements to the ag order in the 
paper. In order to sec explicitly the effects of new operators in the model 111 
2HDM we divide the decay amplitude into two parts. One has the same form 
as that in SM, the other is new. That is, we can write the decay amplitude, 
to the as order, for B — > (f)K in the heavy quark limit as [15, 21, 13] 

A{B 4>K) = < (/)|s7^s|0 >< K\s-i''h\B >, 
v2 

A^A" + A", (37) 
A" = K6C[«3 + al + al-\{a-, + al + <o) + <oa] 

+ye6K* [«3 + a| + - ^(a? + + a^o) + alJ] , (38) 
= -Vt,V;Jar + —[-lKs{ai2 + a[,) + A,, ^ {a,, + a[,)]\39) 

\ nib 2 ITT'b ' 

For the hadronic matrix elements of the vector current, we can write < 
0|s7^6|O >= m^j^el and < K\s-i^h\B >= F^^'' {q^){p%^fj,)^{F^-'^ {q^)- 
Ff'~*^{q^)){m% ~ f^'idff' / (f ■ Here, the coefficients a"'^ in eq.(38) are given 
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by^ 



«3=«3 = C3+- + -— C4F^, 

= C4 + — + [c3(F^ + G^{ss) + G^{sb)) + C2G^{sp) 

b 

+ (c4 + ce) G4>{sf ) + csgG^^s 

f=u 

< = «5 = C5 + ^ + ^^ca(-F,-12), 

a^ = 4 = c,+ | + ^^cs(-F,-12), 

u c Clo agCp 

u c Cg ttsCp „ 



47r TV 



(c8 + cio)- ^fG,p{sf) + cg-{esG^{ss) + ebG^{sb))^^) 

f=u 



where p takes the values u and c, N — 3, Cp — (-/V^ — 1)/2N, and s/ 
G4>{s) = - + -ln-^-G'^(s), 

G^{s)-4:J dx^^{x) J duu{l - u)lii[s - u{l - u){l - x)], 



G^,s,^-2Gl, Gl^ j 

F^ = -121nA_i8 + 4 + 4^, 

/•I 1 — 2x 

fi^ dxg{x)(t)4,{x), g{x)^3- Inx-Sin, 

^ Jo 1 — X 



^The explicit expressions of the coefficients a"''^ have been given in ref. [21]. Because 

there are minor errors in the expressions in the paper and in order to make this paper 
self-contained we reproduce them here, correcting the minor errors. 
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B 



N Ff-^^(0)m| Jo z 
where (f)i{x) are meson wave functions, 

(f)B{x) = Nbx'^{1 - xfExp[- 



r ,,tsM r r dy^, (4i) 

JO z Jo X Jo V 



.2 ^2 



(42) 



with normahzation factor Nb satisfying dx(j)B{x) — 1. Fitting various B 
decay data, lob is determined to be 0.4 GeV [22]. In eq.(42) the asymp- 
totic hmit of the leading-twist distribution amphtudes for and K has been 
assumed. 

The coefficients Oj in eq. (39) are 

„ TDneu 



where 



and 



neu 
U4 



ai2 



ai6 



47r iVc 



C12 + 



C16 + 



15 



TV, ' 



l^' = -121n-^ -6 + C dxgi 
mi, Jo 



X 0MX 



(43) 



(44) 



Tjneu 



X 



rrisX* /4 , rrib 



--(Cn -r -^11 



+(^13 -I- (7(3) Abb 



^,„_._G,(0)UA^iln^-G,(l) 



-21n^G5-GF^(l) 
A* 
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In eq. (45) 



-4(Ci5 + C[,)X^ 



bb 



-8(Ci6 + C[ 



A 



bb 



GFjr 



rribJ V A* > 



GF^s) = dx^^GF(s-ie,x), 
Jo X 

GF(s,x) = / dUnfs-xii], 
Jo 



(45) 



(46) 



with x = 1 — x. In calculations we have set rriu^d = and neglected the terms 
which are proportional to rri^/ml in eq. (45). We have included only the 
leading twist contributions in eq. (43). In eq.(39) a'^ is obtained from by 
substituting the Wilson coefficients Cjs for CjS. In numerical calculations a[ 
is set to be zero because we have neglected C'jS. Wc sec from eq. (45) that 
the new contributions to the decay amplitude can be large if the coupling Abb 
is large due to the large contributions to the hadronic elements of the NHB 
induced operators at the ag order arising from penguin contractions with b 
quark in the loop. 

The decay rate can be obtained [21] 



n2 



2 f2 ttB^K, 



AVfiFt 



m 



2\2 



mlP, 



3/2 



(47) 



where P, 



(1 — m1/m\ — mym%Y — Amfmym%. 
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The time-dependent CP-asymmetry S^k is given by 



a^K{t) = -C^KCOs{AMBot) + S^Ksm{AMBot), 



(48) 
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where 



Here A^^ is defined as 

The ratio {q/p)B is nearly a pure phase. In SM A^^ = e'^^-|-0(A^. As pointed 
out in Introduction, the model III can give a phase to the decay which we 
call (jP^^. Then we have 

^ ^ ^.(2/3+^"^) \A ^S^K = sin(2/3 + (51) 

if the ratio ^| = 1. In general the ratio in the model 111 is not equal to 
one and consequently it has an effect on the value of S^k, as can be seen 
from eq. (49). Thus the presence of the phases in the Yukawa couplings of 
the charged and neutral Higgs bosons can alter the value of -S^k from the 
standard model prediction of St^K — sin 2(3j/^k ~ 0.7. 

5 Numerical analysis 
5.1 Parameters input 

In our numerical calculations we will use the following values for the relevant 
parameters: rub = 4.8 GeV, mc = 1.5 GeV, m* = 175 GeV, A^^^ = 225 MeV, 
2 X 10-^ < Br(B A:,7) < 4.5 x 10-^ d„ < IQ-^^e-cm, = 0.233 GeV, 
/k = 0.158 GeV, fs = 0.18 GeV, Ff-*^{m^) = 0.3. The parameters for 
CKM are = 0.2229, S13 = 0.0036, S23 = 0.0412 and ^13 = 1.02. 
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5.2 Constraints from B Xg^ and neutron electric 
dipole moment (NEDM) 

It is shown in Ref. [14] that the most strictly constraints come from B Xg^ 
and neutron electric dipole moment (NEDM). For completeness, we write the 
formulas as following [24]: 



Br{B^Xs^) r 




Br{B ^ XeePe) 




2 



-|Cr,(m,)^ (52) 

where f{z) = 1 - 8^^ + 8^^ - - 2Az^\nz and Br{h ce-p) = 10.45%. 
The NEDM can be expressed as 



with 



5.3 Numerical results for B Ks(p 

We have scanned the parameter space in model III, in the following we will 
show the results for several specific parameters. 

Figs. 1-4 are devoted to the case in which neutral Higgs boson masses 
are set to be m^jo = 115 GeV, m^o = 120 GeV, ttiho — 160 GeV, which are 
the same with Ref. [19], and consequently Cii(r?iiv-) ^ Cisimw)- Figures 1 
and 2 show the Sfj^x and AS, defined as the S'^x difference with and without 
NHB contributions, as a function of ^ = 6bb + Ou with mH± — 200 GeV. Note 
that there is another allowed region of 9, about —1.2 ~ —0.7, in which S^k 
is about 1. Therefore, we don't present the results in the figures. From the 
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figures we can see that in model III, tlie cliarged and neutral Higgs boson 
contributions can decrease the value of S(j,K down to -0.2, in the allowed 
parameter space. It should be emphasized that the NHB contributions are 
sizable. In fig. 3 and 4, we show the direct CP violation variable C,j,K and 
AC^i^, defined as C<^x difference with and without NHB contributions, as 
a function of 6. It is obvious that C^k can be 8-20 %, i.e., it can be in 
agreement with the data within Icr deviation, while it is zero in the SM. At 
the same time, the NHB contributions can only change C^k by less than 3%. 

Figs. 5-8 [and also in figs 9-11] are plotted for the case in which the 
masses of NHBs have large splitting, m^o = m^^o = 1 TcV 3> m/jO = 115 
GeV, and consequently Cii{mw) is the same order of magnitude, compared 
to Cizijnw)- Now the NHB contributions are as important as those of the 
charged Higgs boson and S^k can reach about —0.6, as expected. 

In order to demonstrate the NHB contributions, in figures 9-11, we show 
S^K as functions of the phases of \hb and As^, and Oss-, and the correlation 
between S^k and C^k-i respectively. It is clear that S^k is sensitive to the 
phases. At the same time, in the range [— tt, tt] of 9},}, and Oss C<j)K changes 
only several percents. There is a strong correlation between S^k and C^k 
and C(i)K is always positive regardless of the sign of jS'^r-, which is opposite 
to that of the central value of measurements. Therefore, if the minus C^^- is 
confirmed in coming experiments the model III 2HDM could be excluded. 
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6 Conclusions and Discussions 



In summary we have calculated the Wilson coefficients at NLO for the oper- 
ators in the SM [except for Qf^ and Qsg], and at LO for the new operators 
which are induced by NHB penguins in the model III 2HDM. Using the 
Wilson coefficients obtained, we have calculated the mixing induced time- 
dependent CP asymmetry S^k, branching ratio and direct CP asymmetry 
CifiK for the decay B (pKg. It is shown that in the reasonable region of pa- 
rameters where the constraints from B — B mixing ,T{h ^ s'j), r[b ^ ctVt)-, 
po, Rb, B — > iJi^n~ , and electric dipole moments (EDMs) of the electron and 
neutron are satisfied, the branching ratio of the decay can reach 10 x 10~^, 
C^K can reach 18% and S^pk can be negative in quite a large region of pa- 
rameters and as low as -0.6 in some regions of parameters. 

Let us separately discuss the two cases: 1) only the charged Higgs con- 
tributions and 2) only the NHB contributions, in addition to the SM ones. 
Without NHB contributions, i.e., in the first case, the charged Higgs contri- 
butions can only decrease S^k to around 0. That is, the model HI can agree 
with the present data, S(f,k = —0.39 ± 0.41, within la error. 

For the second case, our results show that the effects of NHB induced 
operators can be sizable even significant, depending on the characteristic scale 
ji of the process. Due to the large contributions to the hadronic elements of 
the operators at the order arising from penguin contractions with b quark 
in the loop, both the Br and S(j)K are sizable or significantly different from 
those in SM. 

Putting all the contributions together, we conclude that the model HI can 



17 



agree with the present data, S^k = —0.39 ± 0.41, within the la error. Even 
if the S(j,k is measured to a level of —0.4 ±0.1 in the future, the model III can 
still agree with the data at the 2(7 level in quite a large regions of parameters 
and at the lu level in some regions of parameters. As for C<^/<-, our result is 
that it is positive, which is opposite to that of the measured central value. 
Considering the large uncertainties both theoretically and experimentally at 
present, we should not take it seriously. 

Our results show that both the Br and S(f,k (as well as C^k) of S — > (t)Ks 
are sensitive to the characteristic scale of the process, as can be seen 
from eq. (45) and the SM amplitude. The significant scale dependence 
comes mainly from the 0{aa) corrections of hadronic matrix elements of the 
operators Qj, i=ll,...,16 and also from leading order Wilson coefficients Cj, 
i=8g,ll,...,16. However, despite there is the scale dependence, the conclusion 
that the model III can agree with the present data, 5'</,fe = — 0.39±0.41, within 
the lex error can still be drawn definitely. 

Note Added 

We noticed the reference [26] while completing this work. In the ref. [26] the 
mixing induced CP asymmetry S^k in the model III 2HDM is investigated. 
Comparing with the paper, our results on the Wilson coefficients of the op- 
erators which exist in SM at NLO are in agreement. We differ significantly 
from the paper in the neutral Higgs boson contributions included. Further- 
more, we calculate hadronic matrix elements of operators to the Us order 
by BBNS's approach while the paper uses the naive factorization, i.e., at 
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the tree level. Therefore, our numerical results and consequently conclusions 
are significantly different from those in the paper. Even without including 
the NHB contributions our results are also different from theirs due to the 
different precisions of calculating hadronic matrix elements, to which S^pk is 
sensitive. 

When the publication processing the paper, we were aware of Ref. [27] in 
which the LO anomalous dimensions for the mixing of (5ii,i2 onto (53,...,6 and 
Qq are given and those for the mixing of (5i3,...,i6 onto Q7-y,%g given in Ref. [28] 
are confirmed. In this paper these mixings are not taken into account. If in- 
cluding them the numerical results would change but the qualitative features 
of the results would be the same. We shall include them in a forthcoming 
paper on "CP asymmetries m. B ^ rj'Ks and (pKs in a model III 2HDM". 
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Figure 1: 5'<^x as a function of 6* = 9bb + 9tt with /i = 2mb (solid), rrib 
(dashed) andm(,/2 (dotted), where mij± = 200 GeV, \Xu\ ~ 0.03, \Xbb\ = 100, 
^ss — Kc — lOOe"*'^/^. The parameter in neutron EDM expression is 0.03 
[14, 25]. 
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Figure 2: [defined as difference between S^k witli and witliout NHB 
contributions] as a function of 6 with fi = 2mb (solid), (daslied) and mh/2 
(dotted). Other parameters are the same with Fig. 1. 
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Figure 3: C^k as a function of 6. Other parameters and conventions are the 
same with Fig. 1. 
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Figure 4: AC [defined as difference between C^k with and without NHB 
contributions] as a function of 6. Other parameters and conventions are the 
same with Fig. 3. 
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Figure 5: S^k as a function of 9 with ji = 2mb (solid), nib (dashed) and 
mb/2 (dotted), where mH± = 200 GeV, |A«| = 0.03, \\bb\ = 100, A^, = A^c = 
lOOe-*'^/^. Note that the masses of NHB [in Figs. 5-11] are different with 
those in Figs. 1-4. 
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Figure 6: as a function of 9 with = 2mh (solid), rrib (dashed) and mb/2 
(dotted). Other parameters are the same with Fig. 5. 
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Figure 7: C^k as a function of 6. Other parameters and conventions are the 
same with Fig. 5. 
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Figure 8: AC as a function of 9. Other parameters and conventions are the 
same with Fig. 7. 
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Figure 9: S^k as a function of 9bi with fi = 2mb (solid), mb (dashed) and 
mb/2 (dotted), where mj^± = 200 GeV, lA^I = 0.03, \Xbb\ = 100, 9 = 1.15 
and Xcc = A^, = lOOe^'^/l 
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Figure 10: Correlation between Cfpx and Sfpx, other parameters are the same 
with Fig. 9. The outmost two curves correspond to /i = mb/2, the curve in 
kernel is for fj, = mb and the other two curves are for /x = 2mb- 
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Figure 11: S(f,K as a function of 9ss with /i = 2mb (solid), (dashed) 
and mb/2 (dotted), where mj^± = 200 GeV, lA^] = 0.03, |A,,| = 100, 9 = 
= lOOe-''^/^ and A^c = lOOe^'^/l 
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